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Abstract 

The experimental discovery of large — u-r mixing indicates that analogous mixing in 
the charged lepton sector could be substantial. We consider the possibility that if a high 
intensity muon beam, perhaps at the early stages of a muon or neutrino factory, strikes 
a nuclear target, then conversion of some of the muons into tau leptons could occur 
(similar to the conversion of muons to electrons at MECO). Using current experimental 
limits on rare tau decays to bound the size of the relevant operators, we find that a 
50 GeV muon beam, with 10^" muons on target per year, could yield as many as 10^ 
fj, + N —f T + N events per year. Backgrounds could be substantial, and we comment 
on the possibility of detection of this process. 
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In the past decade, the biggest surprise in our understanding of flavor physics has been 
the discovery of large mixing |JJ in the neutrino sector. This large mixing may come from 
diagonalization of the neutrino mass matrix, the charged lepton mass matrix, or both. It is 
quite possible that searches for charged lepton flavor violation will be critical in determining 
the physics of flavor violation. 

The most promising such search is the MECO experiment [2j, in which a high intensity, low 
energy muon beam strikes a nuclear target. The muons are captured, and decay essentially 
at rest. The conversion fi + N ^ e + N will then yield a distinctive 105 MeV electron. The 
experiment promises to achieve the extraordinary sensitivity of a part in 10^'', and in many 
models beyond the standard model, a positive signal would be expected. 

Observations of atmospheric neutrinos indicate that mixing between the muon and tau 
neutrinos is maximal. This gives strong motivation for considering transitions between the 
muon and the tau. Of course, the analogous process to MECO, t + N fJ' + N, is impractical 
due to the short lifetime of the r. However, the inverse process, fi + N ^ t + N might be 
possible. Unlike MECO, this can't occur for muons at rest, but in a higher energy muon 
beam, one can look for such events. Such high energy and high intensity muon beams are 
expected P| at neutrino factories (or early stages of muon factories), in which intensities of 
10^° muons per year and beam energies up to 50 GeV have been proposed. In this note, we 
examine whether the fi + N ^ t -\- N process is feasible at such a neutrino factory. 

The existence of the process fi + N t + N immediately implies that there will be 
muon and tau number violating rare r decays, such as r — /ivr, r — > funr, r fip, etc. The 
non-observation (as yet) of these decays implies an upper bound on fi + N—^r + N. We 
first examine the upper bound on the size of the various operators. 

The relevant operators are of the form [fl T r) (g" F q^) , where F contains various combi- 
nations of Dirac gamma matrices. A detailed analysis of all 48 possible operators, where the 
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g's are any combination of the six quarks and the T consists of (1, 75, 7/^, 7^*75) was carried 
out recently by Black, et al. |3]. They determined the experimental lower bound on A for 
each process, where A is defined by the considering the relevant operator to be 

4-7r 

^(/irr)rrg^). (1) 

For simplicity, we will consider valence quarks only, and will assume that the operators 
are isospin invariant, so that the operators involving w-quarks and d-quarks are the same 
magnitude. Relaxing this assumption will only strengthen our results. Black, et al. find 
that the lower bound on A for F = (1,75,7/^,7/^75) is (2.6, 12, 12, 11) TeV, which come from 
r fin^TT^jT — s> fi7i°,T — > fip and r — > fiir", respectively. Since the bound on the scalar 
operator is the weakest, we will assume that the operator is scalar, and is thus 

%ifir){qq), (2) 

where g is m or and A is greater than 2.6 TeV. No experiment can currently exclude such a 
possibility. When our results are presented, we will briefly comment on the effects of choosing 
one of the other three operators. 

With this operator, we can calculate the cross section for fi + q t + q, and we find that 
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Folding in the parton distribution functions, we plot the cross section for fi + N ^ t + N, 
where is a nucleon, in Figure 1, assuming that the lower bound on A is saturated. For the 
expected beam energy of 50 GeV, the cross section is 0.55 fb. 

With this cross section, we can determine the mean free path. If p is the density of the 
target (in g/cm^), the mean free path is 

A = i ( ^ ) (1.6 X 10^^) meters. (4) 



p \ a 
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Figure 1: The cross section for the scattering fiN ^ t N in units of fb as a function of muon 
energy (GeV) in Lab frame. The sohd (dashed) hne represents the cross section assuming a 
scalar (vector) interaction. 

For a 50 GeV muon beam, there is httle ionization loss over a meter of target, and thus there 
is a probability of approximately 3 x lO^^^p of interacting in a meter of target. With 10^° 
muons on target is a year, this gives 3 x lO^p events per year per meter of target. 

We have assumed that the interaction is scalar. If it is vector, there is a factor of 8 increase 
in the square of the matrix element (in the massless limit), however the lower bound on A 
is 12 TeV instead of 2.6 TeV, leading to a lower event rate. This is also plotted in Figure ^ 
Nonetheless, even here there could be well over 100, 000 events per year. Using pseudoscalar 
or axial vector operators will give similar results. But for the scalar case, and a fairly dense 
target, the event rate could exceed 10'' events per year. 

Although this seems to be a huge event rate, the backgrounds could be severe. Note, 
however that the cross section for tau pair production through Bethe - Heitler production off 
iron nuclei [S] is much smaller than a femtobarn, and the pt distributions are much softer, so 
tau pair production will not be a problematic background. The major difficulty is identifying 
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a clear signature. A typical r energy, for a 50 GeV incident beam energy, will be tens of GeV, 
and thus its decay distance would be a couple of millimeters. One can imagine alternating 
target and scintillator, but many r's will be missed. The only places in which r's have 
been detected are the clean environment of electron-positron colliders, the Tevatron, where 
the signature is large missing transverse energy, and DONUT [6]. The latter used lead and 
emulsions as the target and detection media, and it isn't clear whether the enormous intensity 
of the incident muon beam would blacken the emulsion (this would depend on the beam size 
and whether the emulsion is cycled in and out). This possibility should be investigated. What 
are the specific decay modes that might be observable? The leptonic decays will clearly be 
swamped by backgrounds. The vrz/ decay mode will lead to a monochromatic pion, but unless 
a r track can be observed, the backgrounds for single pions in the intense muon beam will 
also be very large. One could look at rarer decays, such as the three charged pion (or even 
five charged pion) decays, coming at the end of a very short track. Clearly, detection of this 
process will not be easy, but the event rate is high enough that a clever scheme might be 
able to pick out a signal. 

Are there specific models which predict such a large rate for fi + N ^ t + N7 The 
Standard Model, with massive neutrinos, will have mixing between the /i and the r, but this 
mixing is of the order of m^/m^, and is thus negligible. However, there are a wide variety 
of extensions of the Standard Model, including models with very heavy neutrinos, horizontal 
symmetries, left-right symmetry, supersymmetry, extended gauge and Higgs models, etc., 
and many of these do predict such mixing to occur at a substantially higher rate. The effects 
of fi — T mixing can be parametrized by operators of the form of Eq. (1). As noted earlier, 
the biggest rates for fi + N ^ r + N occur if the operator is scalar, as in Eq. (2) (since the 
experimental limits on the operator are weaker), and thus models with flavor- changing scalar 
exchanges are most promising. For example, in R-parity violating supersjTiimetry the 
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superpotential can be written in the form Xij^LiLjE^ + X'^ji^LiQjDk- If the underlying theory 
giving rise to this superpotential gives a hierarchical structure for A, so that Aj23 is large, and 
a non-hierarchical structure for A', then the operator of Eq. (2) can be generated via scalar 
neutrino exchange. If the couplings are of order unity, and the scalar neutrino mass is of the 
order of a TeV, then the operator will be as large as allowed by bounds on r — /ivrvr and the 
rate ioi fi + N r + N will be as large as discussed in the previous paragraph. Alternatively, 
supersymmetric models at large tan/5 can have very large flavor- changing Higgs couplings [8j, 
and that can also lead to similarly large muon to tau conversion. Thus, we see that plausible 
extensions of the Standard Model exist in which r /ivr+Tr" is near its current limit. 

The early stages of a neutrino or muon factory will involve a high intensity muon beam 
with energies up to 50 GeV. In this Brief Report, we have proposed that such a facility 
may be able to substantially improve bounds on — r mixing, or discover such mixing, by 
looking for muon conversion in nuclei to tau leptons. The event rate could be high, although 
backgrounds will be challenging. In view of the large mixing in the neutrino sector, this may 
be a promising place to search for new flavor physics. 

After this work was completed, we became aware of a very interesting paper by Gninenko, 
Kirsanov, Krasnikov and Matveev Pj. They also considered the process fi + N^T + Nata 
neutrino factory. Instead of considering the vertex involving valence quarks, as we did, they 
considered the four-fermi interaction (/2r)('uc), involving production of a charmed quark. This 
has a substantial advantage over our vertex, which is flavor diagonal, because there are no 
experimental constraints on the size of this interaction (since r's can't decay into a charmed 
meson plus a muon). As a result, they had a much higher event rate, and could consider 
the muonic decay of the r. They performed a simulation of the signature at the NOMAD 
detector. What is new in our work? We considered the flavor diagonal vertex, which is more 
tightly constrained by experimental bounds. Our belief is that it is very unlikely for the 
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four-fermi interaction to be purely off-diagonal in the mass eigenstate basis, and thus the 
existence of the vertex considered by Gninenko et al. will generally imply the existence of 
the vertex that we have considered. In that sense, this work is complementary to theirs. 
Clearly, there is sufficient interest in the possibility of mu-tau conversion in nuclei that all 
experimental possibilities should be considered. 

We are very grateful to Jon Urheim for numerous discussions, and we also appreciate 
discussions with Andy Sessler, David Armstrong, Stan Brodsky, Jack Kossler and Jeff Nelson. 
The work of MS was supported by the National Science Foundation under grant PHY- 
0243400, and the work of IT was supported by the Scientific and Technical Research Council 
of Turkey (TUBiTAK) in the framework of NATO - Bl program. 

References 

[1] C. Walter, Nucl. Instr. and Meth. A 503, 110 (2003) and references therein. 
[2] S. W. Molzon, Nucl. Phys. Proc. Supp. Ill, 188 (2002). 

[3] See M. M. Alsharo'a, et al., Phys. Rev. ST Accel. Beams 6, 081001 (2003) for an extensive 
list of references. 

[4] D. Black, T. Han, H.-J. He and M. Sher, Phys. Rev. D66, 053002 (2002). 
[5] V. Ganapathi and J. Smith, Phys. Rev D23, 75 (1981). 
[6] K. Kodama et al., Nucl. Instr. and Meth. A 493, 45 (2002). 

[7] A. Gemintern, S. Bar-Shalom, G. Eilam and F. Krauss, Phys. Rev. D67, 115012 (2003); 
J. P. Saha and A. Kundu, Phys. Rev. D66, 054021 (2002). 

[8] K. S. Babu and C. Kolda, Phys. Rev. Lett. 89, 241802 (2002). 



7 



[9] S. N. Gninenko, M. M. Kirsanov, N. V. Kransnikov and V. A. Matveev, Mod. Phys. 
Lett. A17, 1407 (2002). 



8 



